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T-H , Abstract. Based on time-dependent radiation-hydrodynamics simulations of the evolution of 
' Planetary Nebulae (PNe), we have carried out a systematic parameter study to address the 
non-trivial question of how the diffuse X-ray emission of PNe with closed central cavities is 
^, expected to depend on the evolutionary state of the nebula, the mass of the central star, and 
• the metallicity of stellar wind and circumstellar matter. We have also investigated how the model 
_ predictions depend on the treatment of thermal conduction at the interface between the central 
'hot bubble' and the 'cool' inner nebula, and compare the results with recent X-ray observations. 
Our study includes models whose properties resemble the extreme case of PNe with Wolf-Rayet 
type central stars. Indeed, such models are found to produce the highest X-ray luminosities. 
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1. Introduction 



Over the past decade, roughly a dozen nearby Galactic Planetary Nebulae were found 

■ to be a source of extended diffuse X-ray emission (e.g. Kastner et al. 2008; Kastner 
^ 2009). Images taken from space by the X-ray observatories Chandra and XMM-Newton 

■ reveal that the diffuse X-ray emission is confined to the closed central cavities of round 
and elliptical PNe, indicating the presence of hot gas with temperatures in excess of 

^— I 10^ K (for recent examples see Guerrero et al., this volume). Basically, these findings 
K>" confirm theoretical predictions of wind-heated 'hot bubbles', even though the observed 
X-ray temperatures are much lower (at least a factor of 10) than obtained from simple 
^ spherical models assuming that the wind kinetic energy is thermalized in an adiabatic 

■ shock. However, more elaborate models taking into account the heat flux due to thermal 
' conduction across the central 'hot bubble', show much lower X-ray temperatures, in good 

agreement with observational evidence (Steffen et al. 2008, henceforth Paper V). About 
y—i 75% of the targeted PNe show no evidence of diffuse X-rays, among them old, evolved 
I , nebulae and objects with open structures. Note that a large fraction of PNe also harbor 
^ ' X-ray point sources at the position of their central stars. This type of X-ray emission, 
• ^ r however, is of a different nature than the extended X-ray emission from the shock-heated 
hot wind bubbles which are the subject of the present study. 
^ The number of PNe with detected diffuse X-ray emission is too small to allow the 

■ " " ' derivation of reliable empirical trends, such as relations between stellar wind parameters 
and X-ray properties of the nebula. Likewise, it is difficult to make theoretical predic- 
tions based on simple back-of-the-envelope calculations. In very general terms, the X-ray 
luminosity, Lx, is given by the X-ray emissivity per unit mass and electron density, rj, 
integrated over the emitting volume: Lx = r]{Tc, Z) Mxn^ ^ vi^ci Z) ^x^- The emis- 
sivity depends on the electron temperature. To, and the metallicity of the plasma, Z; the 
mean electron density is proportional to the X-ray emitting mass, Mx, divided by the 
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X-ray emitting volume, ~ Mx Rx^- If we want to know how the X-ray luminosity of an 
object changes in the course of its evolution, we have to know how the product R^^ 
changes with time. This depends mainly on the time evolution of the fast wind of the 
central star: for a constant wind expanding into an environment with a radial density 
profile p oc r^^, it can be shown that Mx oc t, and i?x oc t (Koo & McKee 1992), such 
that the X-ray luminosity decreases as Lx (xt~^. On the other hand, the more realistic 
assumption that wind power and momentum flux increase with time, e.g. M V'^ oc t 
MV (X t^-°'^, leads to an increasing X-ray luminosity, Lx oc t^'^ (Zhekov & Perinotto 
1996). Moreover, properties of the X-ray emission also depend on the thermal structure 
of the hot plasma, which is ruled by the efficiency of electron heat conduction. Clearly, 
detailed numerical models are necessary to make reliable theoretical predictions. 



2. Model calculations 

Using the Potsdam NEBEL code, we model the combined evolution of central star 
(CS) and circumstellar envelope by radiation hydrodynamical simulation in ID spherical 
geometry, accounting for non-equilibrium ionization of the nine most important chemical 
elements. A typical model has a radial extent from 6 x 10'^'* to 3 x 10^* cm (0.0002 to 
1.0 pc). Treated in a physically consistent way, it contains the freely expanding CS wind, 
the inner reverse shock, the 'hot bubble' consisting of the shocked stellar wind, and the 
(double-shell) PN proper, separated from the 'hot bubble' by the contact discontinuity 
and from the surrounding AGB halo by the outer shock. Temporally, the simulations 
start at the tip of the AGB and are advanced well into the white dwarf regime, thus 
covering the formation and complete evolution of the PN. 

Given the radial temperature and density structure of the hydrodynamical model, we 
apply the Chianti code (v6.0.1, Dere et al. 1997, 2009) to compute, for selected instants 
along the evolutionary sequence, the emergent X-ray spectrum, the X-ray luminosity, Lxj 
the characteristic X-ray temperature, Tx, and the surface brightness distribution of the 
X-ray emitting hot central cavity. For details about the model calculations see Paper V 
and references therein. 



3. Simulated X-ray emission: a parameter study 

The work by Steffen et al. (2008) was restricted to models with standard Galactic-disk 
elemental abundances, as given in their Tablet. The present parameter study considers 
the influence of the adopted chemical composition on the resulting X-ray properties. 

Models with scaled metallicities. We have obtained model sequences of different metal- 
licity, Z, by globally scaling the number densities (relative to H) of the 'metals' C, N, 
O, Ne, S, CI, and Ar by a common factor fz with respect to the standard Galactic- 
disk abundance mix, Z ~ fz Zqd . As explained in detail in Schonberner et al. (2010) 



changes of the metal content affect both the cooling properties of the nebular matter 
(higher Z more efficient line cooling), and the wind of the central star (higher Z 
higher wind power). Specifically, we assume here M oc Z"-^^ and Kvind oc Z^-^^, so that 
Lwind = l/2^^Kvind °^ Z'^'^^ . At the Same time, we assume that the evolution of the 
central star, Tcs{t), L^,{t), is unaffected by changes of Z (see Schonberner at al. 2010). 

Our previous studies at standard Galactic-disk composition, based on models including 
heat conduction, had shown that the X-ray luminosity increases steadily with time (Tcff ) 
until the turnaround point (maximum Toff) is reached, after which L*, Lwindj and Lx 
drop rapidly. Essentially, the time evolution of Lx refiects the evolution of Lwind- The so 
far unexplored metallicity dependence of Lx and Tx is shown in Fig.[l] for a central star 
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Figure 1. Evolution of the X-ray luminosity, Lx, as a function of the CS effective temperature 
(left), and of the X-ray temperature, Tx, as a function of the CS wind speed (right), for a central 
star with M» — 0.595 M0. The solid curves connecting data points (symbols) refer to model 
sequences computed with heat conduction HC2 (see Paper V), assuming metallicity scaling 
factors of /z = 3 (triangles), 1 (diamonds), 1/3 (squares), and 1/10 (crosses). Corresponding 
sequences computed without heat conduction are shown by the dashed curves (right panel only), 
omitting the data points for fz = 3, Vwind > 2000 km/s which are numerically uncertain. Stars 
mark the position of observed PNe with normal (H-rich) chemical composition (data M. Guerrero 
2011, priv. comm.). Lx and Tx refer to the X-ray emission between 5 and 40 A (0.3 to 2.5 keV). 

with Af* — 0.595 M0. For given Toff, Lx increases systematically with metallicity (left 
panel). Naively, one might have expected that roughly Lx oc Z^, since both iwind and 
the number of metal ions scale with Z. The numerical results, however, show a weaker 
Z sensitivity, Lx cc ZP,1 < p < 2, where p depends on Tcs- The X-ray temperature 
increases only slightly with metallicity (right panel), in quantitative agreement with the 
scaling suggested by Eq. (14) in Paper V, noting that Lwind « Z. In the absence of heat 
conduction, Tx is significantly higher (Fig.[Tl right panel) and - except for the earliest 
stages of evolution ~ Lx is reduced by roughly 2 orders of magnitude (not shown) . 
We have also investigated sequences of different metallicity for = 0.625 Mq and 
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Figure 2. Left: Evolution of Lx over CS effective temperature, for M* — 0.595 Mq. The two 
thick; curves refer to model sequences with standard chemical composition, the corresponding 
thin curves show the effect of changing the CS wind chemistry from standard (H-rich) to WR- 
type (H-deficient) at approximately identical M{t), Vwind(i), and unchanged H-rich composition 
of the nebula. The upper two tracks (triangles) have M{t) increased by a factor 100 (T4rind(i) 
unchanged) with respect to the standard mass loss rate used with the lower two tracks (squares). 
All sequences have been computed with heat conduction HC2 (see Paper V). The observed po- 
sitions of two PNe with WR-type central stars are also indicated (stars marked BD+30°3639 
and NGC 40; data M. Guerrero 2011, priv. comm.). Right: Synthetic X-ray surface brightness 
profile computed for the uppermost track in the left panel at Tcft ~ 90 000 K. From left to right, 
arrows indicate the position of the inner shock, the chemical discontinuity, and the conduction 
front. Both plots refer to the X-ray emission between 5 and 40 A (0.3 to 2.5 keV). 



0.696 Mq. Qualitatively, the Z dependence of the X-ray properties is very similar to 
the case M* — 0.595 M©. Summarizing the combined Z, M*, t dependence, the models 
predict the highest Lx for metal-rich nebulae with evolved (hot) massive central stars. 

Models with Wolf- Ray et-type central stars. Wolf-Rayet (WR) type central stars are 
characterized by a hydrogen-deficient chemical composition of their atmospheres and 
winds, and by mass loss rates that exceed those of normal central stars by factors 10 to 
100. In a first attempt to model the X-ray emission of PNe with WR-type central stars, 
we have modified our standard simulation (M*=0.595 Mq) in two respects. First, we 
changed the chemical composition of the CS wind from standard (H-rich) to WR-type 



The diffuse X-ray emission of Planetary Nebulae 



105 



(2% H, 42% He, 49% C, 5% O by mass) at the tip of the AGB (i=0). This imphes that 
the 'hot bubble' forms out of H-deficient matter, while the surrounding nebula remains 
H-rich. For consistency, we have updated the NEBEL code with a generalized description 
of electron heat conduction that is valid for arbitrary chemical composition (Sandin et 
al., this volume). The impact of changing the chemistry of the CS wind is shown by the 
difference between the two lower curves in Fig. [5] (left panel). The H-deficient 'hot bubble' 
forms at <=1340y, Tcff =15.4kK, Vwind=260km/s, slightly later than the H-rich bubble 
(1190y, IS.lkK, 165km/s), but reaches much higher X-ray luminosities (up to a factor 
100) during the early evolution (low Tes). Obviously, the X-ray emissivity (per unit mass) 
is much higher for WR-type matter. Later, the difference in Lx tends to vanish, since 
the bulk of the hot gas is increasingly made up of evaporated H-rich matter (see below). 
As a second step, we have increased the CS mass loss rate globally by a factor 100. This 
delays the formation of the 'hot bubble' both for normal (1620 y, 19.3 kK, 325km/s), and 
WR-composition CS winds (2160 y, 25.4 kK, 720km/s). As seen in Fig.[2] (left, upper two 
curves), such models easily exceed the X-ray luminosity measured for BD+3G°3639, a 
well studied FN with a [WC]-type central star. Again, corresponding models calculated 
without heat conduction typically have lower Lx by roughly a factor 100, except at the 
very early stages where heat conduction is negligible. Finally, Fig. [2] (right) shows the 
X-ray surface brightness distribution of the test model with WR-type wind composition 
and scaled- up mass loss rate at a late stage of evolution (Tcff w 90 000 K). The 'hot 
bubble' is divided into two chemically distinct parts; the inner region (r ^ 4 x lO^'' cm) 
consists of H-deficient matter, the outer shell consists of H-rich matter that evaporated 
due to heat conduction. In X-rays, the H-deficient core of the central cavity shines much 
brighter than the surrounding H-rich layer (cf. Jacob et al., this volume). 

4. Conclusions 

The observed X-ray luminosities of PNe with normal H-rich composition can be ex- 
plained by hydrodynamical models with heat conduction, either assuming a standard 
Galactic-disk metallicity, .^gd, and a range of central star masses between Af, = 0.565 
and 0.7 Mq, or by a fixed M^, = 0.625 Mq and a range of metallicity scaling factors 
between fz = 1/3 and fz = 'i- The highest X-ray fluxes are expected to originate from 
PNe with WR-type central stars, due to their powerful and H-deficient winds. Models 
with suppressed thermal conduction (magnetic fields) may explain X-ray non-detections. 
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